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Section 1 


Introduction 


This report presents a scalar/vector potential formulation 
for viscous compressible unsteady flows around complex geometries 
such as complete aircraft configurations. 

Several scalar/vector potential methods are available in the 
literature for incompressible flows: these methods are reviewed 
in this Section (with particular emphasis on that proposed by 
this author). The extension of scalar/ vector potential method to 
compressible flows is presented in Section 4. For the sake of 
clarity and completeness some classical results on fluid and 
thermodynamics are presented in Section 2 and 3. It should be 
emphasized that shock waves and turbulence are not addressed in 
this report. 


1.1. Scalar/Vector Potent ia 1 M ethods for Incompress ible Viscous 

Flows 


A review of the state of the art of scalar/ vector potential 
methods for viscous flows is presented here. (This author is not 
aware of any method for compressible flows, and thus this review 
is limited to incompressible flows.) In order to discuss the 
advantages of the approach, the primitive variable approach is 
also briefly reviewed. Incompressible turbulent viscous flows 
are governed by the continuity and Navier-Stokes equations. 
Excellent reviews of the state-of-the-art are given for instance 
in Refs. 1 and 2 and are summarized here. For the sake of 
conciseness, boundary layer formulations are not included in this 
review. 

There are two basic approaches to the solution of these 
equations: solution in the primitive or physical variables (i.e., 
velocity and pressure) and solution in the scalar-vector 
potential variables (also called the vorticity—potential method, 
a three dimensional extension of the two-dimensional vorticity 
s tr e am— f unc t ion method). The relative advantages of the two 
approaches are considered here. 

In the primitive-variable formulation there are three basic 
approaches to the numerical solution of the problem: finite- 
difference, finite-element and nodal methods (see Refs. 1-3). 
The major advantage of working with the primitive variables is 
the simplicity of the equations and the fact that the unknowns 
have physical meaning. However such formulation has considerable 
drawbacks. The use of the Navier-Stokes equations requires the 
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solution in the whole physical space, which may be prohibitive in 
terms of computer costs. Another difficulty connected with the 
primitive-variable formulation for incompres ssible viscous flows 
is the lack of evolution equation for the pressure: the method of 
the artificial compressibility introduced to remedy this problem 
is not fully satisfactory (see Ref. 1). However such a problem 
does not exist for compressible problems, which are the main 
objective of this report. 

Next consider the scalar/ vector potential approach. In two 
dimensions the advantages of the vorticity/stream-f unction method 
over the primitive variable approach are well known. One 
advantage is that the continuity equation is automatically 
satisfied. However, a more important advantage, often ignored in 
literature is that the solution for the equation for vorticity is 
limited to the computional region of the boundary- layer/ wake for 
attached flow, while the equation for the stream function is a 
Poisson's equation which can be transformed into an integral 
equation (also limited to the rotational region). The 
implications is that the exact solution of the Navier-St oke s 
equations can be obtained by studying only the rotational region 
(i.e., boundary layer and wake for attached flows or boundary 
layer plus separated flow). In other words, the 
vort icity/stream-function method eliminates all the disadvantages 
of the formulation in the primitive variables. 

To the contrary of a commonly held idea, the 
vorticity/stream-f unction approach may be extended to three- 
dimensional flows. Such a generalization is referred to as the 
vort icity/potent ial method or the scalar/vector potential method. 
The velocity vector v is given by the general theorem 

v = grad $ + curl A 1.1.1 

where cp is a scalar potential (harmonic for incompressible flows) 
and A is a vector potential. The method is classical: although 
Lighthill (Refs. 4 and 5) is the standard reference for this 
approach, according to Lamb (Ref. 6), who gives a theoretical 
outline of the formulation, this concept was first introduced by 
Stokes (Ref. 7) and later refined by Helmholtz (Ref. 8). The 
decomposition is not unique and depends upon the boundary 
conditions on A. This yields the possibility of different 
'versions' of the same basic methodology: this issue has been 
examined very carefully in Ref. 9, which includes an excellent 
review of the theoretical works on this issue. The theoretical 
foundation for their work is to be found in the work by Smirnov 
(Ref. 10). Important in this respect is also the works of 
Hiraski and Heliums (Refs. 11 and 12), and Ladyzhenskaya (Ref. 
13) . 


The formulation involves solving for the velocity in terms 
of the vorticity using the law of Biot and Savart. The vorticity 
is obtained by analyzing the vorticity transport equation, and 
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determining the vorticity produced at the surface due to the no- 
slip condition. Since the advent of the computer era, this 
formulation has been used by several investigators (Refs. 14-24) 
including the author (Refs. 24-26). Particularly good results 
have been obtained by V/u and his collaborators (Refs. 13, 21 and 
23) . 


1.2. A Scalar/Vector Potential Method for Incompressible Flows 

As mentioned above, all the scalar/vector potential methods 
are based on the classical decomposition of a vector field into 
an irrotational component and a rotational (solenoidal) one. The 
highlights of this version of the method introduced by this 
author and his collaborators in Refs. 25 and 26 are presented 
here. 


The problem is governed by the averaged Navier-Stoke s 
equations for incompressible flows, 

Dv 1 1 * 

— = grad p + - p V v 1.2.1 

Dt p p 

and the continuity equation 


div v = 0 1.2.2 

The boundary conditions are 

v = v B 1.2.3 

where Vg is the velocity of a point on the surface of the body 
and (for a frame of reference fixed with the undisturbed fluid) 

v = 0 1.2.4 


at infinity. 

The method is based on the classical decomposition theorem 

v = grad <f) + curl A 1.2.5 

where <f> is a scalar potential and A is a vector potential, 
related to the vorticity £ by the relationship 

V*A = 1.2.6 

The vorticity is given by the third vortex theorem (obtained 
by taking the curl of Eq. 1.2,1) 

D Z - ix 2 - 

— = (£• grad) v + - V ? 1.2.7 

Dt p 

The equation for the potential, obtained by combining Eqs. 
1.2.2 and 1.2.5, is 
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1 . 2.8 


v* 4* = 0 

Equations 1.2.5 to 1.2.7 are fully equivalent to Eqs. 1.2.1 and 
1.2.2 and are much easier to solve. First Eq. 1.2.6 yields 

A= /// y 4^ dV 1,2,9 

whereas Eq. 1.2.8 can be solved using integral equation methods 
(also known as panel methods). The numerical formulation is 
given in Ref. 26 . 

Finally a brief assessment of the proposed method is given 
in the following. The advantage of the vort ic i ty/po tent ial 
method over the primitive-variable approach has already been 
discussed above (elimination of problem due to lack of pressure- 
evolution equation, etc.). Here it is important to emphasize 
again that (1) the method is fully equivalent to the solution of 
the Navier-Stokes equation and (2) that the solution of Eqs. 
1.2.7 and 1.2.9 requires a (finite difference or finite element) 
grid limited to the nonzero-vorticity region (i.e., boundary 
layer and wake region for attached flows). However, the main 
advantage of the formulation presented above is the fact that it 
may be extended to three-dimensional flows. Such an extension is 
presented in Section 4 of this report. As mentioned above, this 
is believed to be the first time that such an extension has been 
presented. 
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Section 2 

Foundation of Continuum Mechanics 


For the sake of clarity and completeness the derivation of 
the equation of continuum mechanics from fundamental principles 
is outlined in this section. This derivation is classical and is 
similar to the one presented in Serrin7 Application to fluids is 
given in Section 3. 


2.1, Bas ic Pef init ions 

Consider a Cartesian frame of reference in a three- 
dimensional Euclidean space. Let: 


x=f (£ a , t) 


2 . 1.1 


be the (vector) function relating the Cartesian coordinates x v , 
(at time t) of a material point identified by convected 
curvilinear coordinates The coordinates £ a could for example 

coincide with the values of x at an arbitrary initial time t=0. 
The function f is assumed to have an inverse 


£ a = F a (x, t ) 


2 . 1.2 


so that there exists a one— to— one correspondence between the 
point x at time t and its curvilinear coordinates. It is assumed 
that the flow is smooth: in particular surfaces of sharp 
variations in the velocity, i.e., wakes and shock waves, are not 
included in the formulation (see SerrinV w* 226-228). 


An arbitrary quantity g, function of x and t, is also a 
function of and t. The following symbols will be used 

= li. \ 

2t I x - const 


2.1.3 


D g ± 3g 
Dt dt I = const 


I 


2.1.4 


Dg/Dt is called the material (or substantial) derivative of g. 
Note, that by using the chain rule: 

Bl = d A +1 2i| 

Dt dt i dx. dt L* „ . f 

* S - Const 

or 

d 8 _ d S . - , 

Vt dt 8 8 

where v is the velocity of the material point 

i 3x‘| OX 


2.1.5 


2 . 1.6 


v = 


Dt 


= const 


2 . 1. 7 
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2 . 2 . Fundamental Conservat ion Pr inc inle s 


The motion of the continuum is assumed to be governed by the 
following fundamental principles: 


Conservation of Mass 


m. 


pdV = 0 


Conservation of Momentnm 


JFJJJ/ Vdv ' 


/£ pJdv * § tdr 


Conservat ion of Angular Momentum 


M. 


pxKv dV 


■ III pixfdv + £ 


2 . 2.1 


2 . 2.2 


2.2.3 


p«XfdV + ©*xt d^ 

Conserv ation of Energy (or the first law of thermodynamics. 
Ref. 2, p. 177) 


illl 


p(e + jv.v)dV 


“/// Pf-vdV + 

6 


Qdf 2.2.4 


In Equations 2.2.1, 2.2.2, 2.2.3 and 2.2.4 the volume V„is 
an arbitrary material volume (i.e., by definition a volume which 
moves with the continuum particles), v is the velocity of the 
fluid particles defined by Eq. 2.1.7, p is the density (mass per 
unit volume), f is the force per unit mass acting on the fluid at 
a point of the volume V^, t is the stress vector (or force per 
unit surface area) acting on the continuum at a point of the 
surface «5 surrounding V H and h is the heat flux supplied by the 
volume Vj^. 

It may be worth noting that, in postulating the above 
principles as the governing equations of the motion of the 
continuum, it is implicity assumed that we are dealing with 
simple nonreacting species; multiphase flows, and/or chemical 
reactions are not included in the formulation. Also 
electromagnetic phenomena are not considered in this report. 


A considerable amount of results can be obtained as a 
consequence of the above set of principles. The rest of this 
section is devoted to the derivation of such results. 


2.3. J acob ian of Transformation 

Note that, for any arbitrary funtion g, 

JJ/ V gdv ■ /// ’* j d?,d?,d *' 

V r1 V 0 

where V 0 is the image of in the 5* space, J is the 
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Jacobian of the transformation 

j = det (— ) > O 

For future reference, note that: 

— = J div v 
Dt 


2.3.2 


2.3.3 


In order to prove this, note that if X? is the cofactor of 
dr'/ 3£ a so that 

•7 3x1 Y «. 

L A = J6 V 2.3.4 

3 §“ J i 

then, 

2£.2£. ( 2si>x‘-2 2s.V 

Dt Dt ; ;< jj* i 

= 2j X. = Z J - J div V 2.3,5 

?jV J ; 2)X‘ 


2.4. Time Derivative of Integrals 


Note that if V^j is a volume moving with the fluid, then the 
boundaries of its image V 0 in the £ a space are time 
independent. Therefore, using Eqs* 2.3.1 and 2.3.3 one obtains: 


£!•«- J/i'if 


+ gdiv v)dV 


2.4.1 


77 gdV = — (gDd^d^d^ 

Cifc JJJ J JJ Dt 




= jjj + gdiv vJJd^d^dS* 2.4.2 


Also, note that 


+ gdiv v = -h <r) 


hence, applying the divergence theorem 


Jjj div wdV = <jj> w^nde' 

(where n is the outer normal to the surface c) one obtains 


dt Jff gdV= jfl^ v + f § "" j6 


2.4.3 


2.4.4 


2.4.5 
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which states tnat the rate of change of the volume integral is 
due so the rate of change of the integrand over a fixed volume 
pins the flux over the boundary surface. 


2_i5j Continuity Equation 

Consider the principle of conservation of mass, Eq. 2.2.1. 
Using Eq. 2.4.2, Eq. 2.2.1 yields (noting that Vjj is an arbitrary 
volume) 

Dp 

— + pdiv v = 0 2.5.1 


or, according to Eq. 2.4.3, 


dp 

at 


+ div 


(pv) 


= 0 


2.5.2 


Also combining Eq. 2.3.3 and 2.5.1 



0 


2.5.3 


which can oe obtained directly trom Eq. 2.2.1, using Eq. 2.3.1 
and 2.4.2 and noting tnat the volume Vo is arbitrary. Equation 
2.5.3 yields: 

pJ = constant = p 0 J 0 (following particle) 2.5.4 


Equations 2.5.1, 2.5.2, 2.5.3, and 2.5.4 are four different forms 
of the continuity equation. Note that, using Eq. 2.5.3 one 
obtains, for any function g, 

£IJJ 4v 2 - s - s 

For, V* 

£ If! p,dV ’JJI §f. ? 1 2 - 5 - 6 

V*. >/ 

n o 

2 .6 . Stress Tensor 


Consider the principle of conservation of momentum, Eq. 
2.2.2. applied to an infinitesimal volume V rt . Let Z be a typical 
length of the volume V so that 

M 

V - 0( Z * ) 2.6.1 

n 

whereas 

6 = Oil*) 2 . 6.2 
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Letting C go to zero, one obtains, in the limit 


t 


tde = 0 


2.6.3 


Assume that Vjj coincides with an infinitesimal Cauchy's 
tetrahedron, i.e., a tetrahedron with the origin at an arbitrary 
point x and three faces parallel to_the coordinate planes (i.e., 
having outward unit normals "* 1 * an< * *^e ^ ace 

in the first octant with normal n (see Figure 2.1). In this case 
Eq. 2.6.3 yields 

t(*)d6 + + t(-Ii)d<r* + t(-I»)d6, = 0 2.6.4 


or setting t(-i^) = -t(i k ) = -tj. and noting that doj. - n^de . 
(where n k are the components of the unit normal n) one obtains 

t = tj.ni t t*n* + t »n j 2.6.5 


i . e . , 



where Tj^ is the kth component of tj=t(ij). 


2 . 6.6 


The above result can be stated as follows: the forces acting 
on three coordinate surfaces through a given point defines a 
tensor (stress tensor), with components the force acting on 
any surface normal to a given direction n is dependent upon these 
quantities through Eq. 2.6.6, which, in tensor notations, may be 
rewritten as 


t = n.T 

2.7. Cauchy* s Eauat ion of Motion 


2.6.7 


Note that, according to Eq. 2.6.6, and using the divergence 
theorem 


6“ * 

-2 




3Tjfcr 

u X ; 


■ fff « 


div T dV 


where div T is a vector defined as 


div T 

j r k 


2.7.1 


2.7.2 


Combining the conservation of momentum, Eq. 2.2.2 with Eq. 2.7.1 
and using Eq. 2.5.5 one obtains 

dt 


pvdV 


V 


= 111 P ot <lv = Ilf (p * + div * )dV 


2.7.3 






M 


2.5 




Figure 2.1 Cauchy tetrahedron 


2.6 



2.7.4 


or, noting that V M is arbitrary, 

Dv 

p— = pf + div T 

which is called Cauchy's equation of motion (Serrin, 7 p. 135) or 
dynamic equilibrium equation. 


2.8. Symmetry of Stress Tensor 

Note, that according to Eq 2.6.6, 

#*«*”■? 2JJJ £ ( ” T jk I k ldv 

(5 ' g *j' t 

[((* , £: T jkik' IV *£ 


‘//I 


„ <*.• 

V M 

x xdiv T dV + 




ijXik T jk dV 


\J 


iv dv 


2 . 8.1 


v n 

Combining the conservation of angular momentum Eq. 2.2.3, 
with Eqs 2.5.5 and 2.8.1 one obtains (noting that vxDr/Dt = vxv = 
0) 

jr Iff p 5 ” dv ■ Iff (£<”v)dv = /// p x „£i dv 

^ V M V M 

= IJJ p3xf dV + IjJ x xdiv T dV + T fff T. k T.xl k d/ 2.8.2 

V* V ( V \ 

or, using the equilibrium equation, Eq. 2.7 .4, and noting that V 

• * t . • j i 


is arbitrary 




T jk V*k = 0 

J ' ** 


or 


T jk - T kj 


2.8.3 

2.8.4 


which shows that the stress tensor is symmetric. 


2 _JL. Energy Transfer Equation and Virtual Work Principle 

Consider the equilibrium equation, Eq. 2.7.4. Taking the dot 
product of both sides of the equation with v and integrating over 
an arbitrary volume V one obtains: 


Note that, using Eqs. 2.6.6, 2.7.2, and 2.8.4, 


II{ ph ™ * III 

v. 


div T-vdV 


2.9.1 


i 




2.7 



JJj div T.v dV = Z jjl 


ILl k t r 


M 




= 2 


2 fJl 


3Tix 

c) *.• 


i k v dV 


dV = 


M 


/f( - 1 llhjijfj dV 

-J | dff IT vfr + ^ dV 

* \ / J 


-i 


t • v d£ 


- //( ™ 

Y,V\ 


dV 


where D is a symmetric tensor with component 

- j<— * li> 

J 2 C>x 3 k l 


2.9.2 


2.9.3 


^7 


J ^ 

and is called the deformation tensor (Serrin, p. 139) or strain- 
rate tensor (see also section 2.B). 


The last integral in Eq. 2.9.2 is called dissipation ter m 
(Ref. 27 , p. 138), and is the work per unit time done by the 
internal stresses. 


Combining Eqs. 2.9.1 and 2.9.2 
volume : 


yields if V is a material 
f M 




6 


t.v de 


2.9.4 


which states that the time derivative of the kinetic energy of a 
material volume is equal to the work per unit time of the volume 
and surface forces diminished by dissipation term. 


The above result is called energy transfer equation or 
mechanical energy equation and relates 'mechanical energy terms' 
to the dissipation term. 

It may be worth noting that if Eq. 2.7.4 multiplied 
(internal product) by 65 instead of v one obtains the virtual- 


work equation 
principle ) 

(which here is not 

assumed as an independent 

/// dv ■ /// p7 ' 6 ° dv * # 

t • So do -jjj T:SD dV 2.9.5 


v M cr 

v 

where 


> A 


SUjv — — 5ut. 2.9.6 

dx. * 


J 
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2.10. Thermodynamic Energy Equation 


Comparing the principle of conservation of energy, Eq. 2.2.4 
with the energy transfer equation, Eq. 2.9.4, one obtains: 


M. Ill ^ = III P ^ V ' // ™ dV - 

v* v* c 


2.10.1 


which relates thermodynamic energy terms to the dissipation term. 
For this reason Eq. 2.10.1 is here called the integral 
thermodynamic energy equation. 


It nay be worth noting that Eq. 2.10.1 is often referred to 
as the first principle of thermodynamics. Here we assune Eq. 
2.2.4 (conservation of energy, as a fundanental principle. Eq. 
2.10.1 is a consequence of Eq. 2.2.4, not a fundanental 
princ iple . 


2.11. Heat Flux Vector 

Note that in Eq. 2.10.1 if £ is a typical size of the volume 
V and t goes to zero then the surface integral is of order 
whereas the other terms are of order l*. Hence if V M is an 
infinitesimal Cauchy's tetrahedron (see Figure 2.1) one obtains 
(see Section 2.6) that: 

Q(n) =Q(ii)ni + Q(i a )n* +Q(ij)nj 2.11.1 


or 

Q = q-n 2.11.2 

where q is a vector with components: 

q k =Q ( i k ) 2.11.3 

equal to the heat flux per unit area through the surface normal 
to I t . 

2.12, Thermodynamic Energy Equation in Differential Form 

Combining Eq. 2.11.3 with the thermodynamic energy equation 
Eq. 2.10.1, applying divergence theorem and noting that the 
volume is arbitrary yields: 

P Dt = ™ “ diV * 2.12.1 

which will be called here the differential form of the 
thermodynamic energy equation. 'According to Truesdell this 
equation should be attributed to C. Newman* (SerrinV p.177) 
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2 . A . A Convenient Expres sion for the Ac celera t ion 
Consider 


Note that 


and set 


so that 


Dvfc 3v k V ^ v k 

= +Z*v* 

Dt 3t j J 3xj 




3x k 


Dv 

Dt 


dv 1 a - - 

— + -grad v - vxC 


and Cauchy's law of motion may be written as 

1 

2 


dv - 2 

— + -grad V 

at 


vx£ = f + idiv T 


2.A.1 


2. A. 2 
2. A. 3 

2. A. 4 


2. A. 5 


2 .B. Kinematics of Deformation 

The relationship between the vorticity and the angular 
velocity of a fluid element surrounding a point x 0 is introduced 
in this section. 

Consider she Taylor series of the velocity field 
arbitrary fixed point x 0 in the flow field 

vj(x) = ViUo) + 2. — ^-rj. + 0(r a ) 

3 J k 

where 

r = x - x 0 2 .B.2 

Equation 2.B. 1. may be rewritten as 

v = vo + r . Tl + r . D 2.B.3 

where the tensor gradv has been decomposed into its symmetric and 
antisymmetric parts 


about any 
2.B.1 


gradv = D 


2.B.4 


where D is the (symmetric) deformation tensor with components 

i 3vj 8v k 

D ik = + T~ ) 2.B.5 

J 2 ox k dxj 
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(see Eq. 2.9.3), whereas Ji. is the (antisymmetric) 
tensor with components 


n 


i dvi 


kj 


3v k 


) 


rotation 


2.B.6 


Note that 

where 

Also 


where 


_/7 k j r k ij = uxr 

® = (rt*» , 12) = ^curl v = Lz, 

r.° = 2 D kj r k ij 
J 


= -grad D 


D = ^ r j D jk r k = 


2.B.7 

2.B.8 


2.B.9 


2.B.10 


is called strain rate qnadric. Hence Eq. 2.B.3 may be rewritten 
as 

v = v 0 + Soxr + vjj 2.B.11 


with 


1 

V D = j grad D 


2.B.12 


which indicates that the motion of a fluid element aronnd a point 
x 0 can be decomposed in translation (with velocity vo = v(xo)), 
rotation (with angular velocity w 0 equal to half the vorticity, 
curl v, at x 0 ) and deformation (with deformation velocity equal 
to grad D/2). 
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Section 3 


Entropy and the Second Law of Thermodynamic s 


The fundamental equations of continuum mechanics were 
derived in Section 2, from the fundamental principles of 
conservation of mass, momentum, angular momentum and energy in 
the form of Eqs. 2.2.1, 2.2.2, 2.2.3, and 2.2.4. In this section 
the formulation is carried further by assuming that the continuum 
be a fluid: the equation of state for a fluid is postulated along 
with the second law of thermodynamics. The formulation is again 
classical and the presentation given here is quite similar to the 
one of Serrin’( pp. 172-178 and 230-241). 


3.1. Starting Eouations 

Starting from the basic principles of conservation the 
following equations were derived in Section 2: 


Continuity equation (Eq. 2.5.1) 

+ p div v = 0 
Dt 

Cauchy's equation of motion (Eq. 2.7.4.) 

Dv - — 

p — = pf + div T 
Dt 

Symmetry of stress tensor (Eq. 2.8.4) 


3.1.1 


3.1.2 



Thermodynamic energy equation (Eq. 2.12.1) 

De 

p— — = T:D - div q 


3.1.3 

3.1.4 


where the deformation tensor, D, is defined by (Eq. 2.9.3) 

j dvt 0Vi 

V - 2*51? + 4’ 

In Eq. 3.1.2, the _force per unit mass f is prescribed. It is 
apparent that if T and q were known, then Eqs. 3.1.1, 3.1.2 and 
3.1.4, with appropriate initial conditions and boundary 
conditions, could be used to obtain p, v and e. In other words, 
the above equations cannot be used to solve the problem unless 
suitable constitutive relations (relating T and q to other 
quanitities) are available. The objective of this section is to 
show that the second law of thermodynamics (introduced in Section 
3.4) puts a constraint (Eq. 3.4.4) on the general nature of the 
constitutive equations. The actual constitutive equations are 
introduced in Section 3.B. 


3.1 



3.2. Equation of State 


Our continuum is assumed to be a fluid, i.e.» a single-phase 
system which is described by two state variables, for instance 
internal energy e and density p introduced in Section 2. All 
other variables are assumed to be functions of the first two. 
Any two variables however may be chosen as the main variables: 'A 
particularly elegant formulation of this relation is that of 
Gibbs' (SerrinJ p.172), which is followed here: the two main 

state variables are chosen to be the entropy (which is assumed 
as a primitive concept, like the energy), and the specific volume, 
r = l/p. The fundamental state equation is some definite 

relationship giving the internal energy as a function of entropy 
and specific volume, of the type 

e = e(S, r) 3.2.1 


In addition, the (thermodynamic) pressure and the temperature 
are defined by 


P = 
T = 


3e 

0T 

0e 

as 


3.2.2 

3.2.3 


(it is assumed, of course, that p and T are greater than zero). 
Since e and p=l/r have already been introduced, Eq. 3.2.1 could 
be thought of as the (implicit) definition of entropy. Note also 
that the pressure is introduced as a thermodynamic quantity (see 
also Eq. 3.3.4) rather than a mechanical one (i.e., force per 
unit surface). 


The differential of the internal energy is 

de = TdS - pdr 

and, accordingly, the material time derivative is given by 

De _ DS Dr 

Dt Dt 1 Dt 


3.2.4 


3.2.5 


3.3. Entropy Evolution Equation 


Combining the thermodynamic energy equation, Eq. 3.1.4, 
with Eq. 3.2.5 one obtains 


DS Dr 

pT pp — = T:D - div q 

Dt Dt 

Note that r=l/p and therefore, according to Eq. 3.1.1 

Dr 1 Dp 

p— — = div v 

Dt p Dt 


3.3.1 


3.3.2 


3.2 



3.3.3 


Combining Eq. 3.3.1 and Eq. 3.3.2 yields 

DS _ - _ 

PT— = pdiv v + T:D - div q 

Equation 3.3.2 may be rewritten in simpler form if we introduce 
the tensor, V, defined by 


V = T + pi 3 .3 .4 

where p is the thermodytwmic pressure (Eq. 3.2.2) and I is the 
unit tensor. Using Eq. 3.3.4, Eq. 3.3.3 may be written as 

pT^ = - div q 3.3.5 

where 


= V:D = T:5 + pdiv v 


3.3.6 


Equation 3.3.5 will be referred to as the entropy- evolution 
equation. It should be noted that Eq. 3.3.5 was obtained from 
the total energy equation (Eq. 3.1.4), by replacing the internal 
energy, e, with its expression obtained from the fundamental 
state equation. Therefore Eq. 3.3.5 is fully equivalent to the 
principle of conservation of energy. 


3.4. Second Law o f Thermodynamic s 

The main postulate introduced in this section is the second 
law of thermodynamics: 

77 (j[ P Sdv > ~ jf ^ q-n dir 3.4.1 

at J v n e 1 

where \^is a fluid— volume (i.e., by definition, a volume moving 
with the fluid particles, see Section 2.2). 


Equation 3.4.1, may be rewritten in a much more interesting 
form, by noting that, using the entropy-evolution equation, Eq. 
3.3.5 (see also Eq. 2.5.5), one obtains: 




III ' f ‘ i dlv 5,dV 

v n 


Hence the second 
rewritten as: 



<£ 1 . 

q-grad T)dV 

T T > 

of thermodynami 



cs, Eq. 3 


d 6 


.4.1, 


3.4.2 


may be 


3.3 



d_ 

dt 


fff pSdV +|lq.5 de = 
V„ * 


■I 


<1 1 

( -q.grad T)dV > 0 

T T * 


3.4.3 


or. 


since V is arbitrary. 


<%> “ — < q.grad T > 0 


3.4.4 


(this equation will be used to discuss the constitutive 
relations. Section 4.A). 


It is instructive to rewrite the total energy equation, Eq. 
3.1.4, as (see Eq. 3.3.6) 

De 

p — + pdiv v = - div q + <f> 3.4.5 

which shows on the right hand-side of the dissipative terms, 
whereas the terms on the left hand 1 side are nondissipat ive. Both 
terms are equal to pT DS/Dt (Eq. 3.3.5), which indicates that the 
entropy is the 'link' between dissipative and nondissipat ive 
terms. 


3 .5 . Thermodynamic Pressure and Mechanical Pressure 

It is worth noting that Eq. 3.3.4 may be rewritten as 

T = - pi + V 3.5.1 

where p is still the thermodynamic quantity originally defined by 
Eq. 3.2.2. However it is apparent that in Eq. 3.5.1, p assumes 
the role of mechanical pressure (force per unit area). This is 
quite clearly indicated by the fact that, in Eq. 3.4.4, (see 
Eq. 3.3.6) is affected by V, not by T. Therefore the second law 
of thermodynamics itself suggests that T be given in the form of 
Eq. 3.5.1 (with V responsible for the dissipative effects), i.e., 
that the nondissipative part of T be given by -pi, where p J.s the 
thermodynamic pressure (in particular for perfect fluids, T==— pi). 
This point, clarified further in Section 4£, is somewhat obscure 
in the literature where Eq. 3.5.1 is introduced as an independent 
assumption. 


3.6. Enthalpy 

Note that, using Eq. 3.5.1, Cauchy's equation of motion may 
be rewritten as 

Dv - 1 1 — 

— = f grad p + -div V 

Dt p P 


3.4 


3.6.1 



In discussing Bernoullian theorems, it is convenient to express 
(1/p) grad p in terms of a conservative (i.e., exact 
differential) and nonconservative part: note that according to 
Eq. 3 .2.4 

1 , . 

—dp = rdp = d(rp) - pdt 
P 

= d(e + rp) - TdS 3.6.2 

and hence 

—grad p = grad h - Tgrad S 3.6.3 

P 

where the enthalpy h is defined as 

h = e + — 3 .6 .4 

P 

Accordingly Eq. 3.6.1 may be rewritten as 

Dv 1 _ 

— = f - grad h + Tgrad S + -j^-div V 3.6.5 

Note that for barotropic flows 0 

3.6.6 

Here Eq. 3.6.3 has the same 'role' that Eq. 3.6.6 has for 
barotropic-f low formulations (for isentropic flows Eq. 3.6.3 and 
Eq. 3.6.6 coincide since, in this case dh=dp/p). 


1 r dp 

—grad p = grad / — 
P i P 


3 .A Dependence of Internal Energy on Specific Volume 

A classical result on 3u/3x for T = constant is derived here 
for use in Section 4.A. Equation 3.2.4 may be rewritten as 


dS=±de +|dr 


or, using T and t as fundamental variables 

1 3e 1 3e p 

dS = — — dT + — — dr + — d i 
T 3T T 3t T 


Hence, 


3S 1 3e 

3T -- T 3T 


3.A.1 


3.A.2 


3.A.3 


as 

3r 


13e 

= -(— + p) 

T 3t 


3.A.4 


where it is understood that the partial differention with respect 
to T is performed with x = const (and viceversa).Theref ore 


3.5 



3. A. 5 


or 


Equation 


whenever, 


a*s _ a i ae i a%_ 

9t 0T “ at 1 ! ar ” T 0T0T 

- 8 * s _ £_ ( i £® P » 

= 0T0r = 0T T 0 t + T ' 

1 0e . i > 0*e 0p 
T a 0 T T^dTdv 0T 


ii| - t!£| - p 

8rl x 8 TI t 
3. A. 6 indicates that 



for r = constant, p is a linear function of T. 


3.A.6 


3.A.7 


3.6 



Section 4 

Vector and Scalar Potentials 


The fundamental equations governing the motion of a fluid 
were derived in Sections 2 and 3. In this section the 
fundamental equations will be rewritten to obtain a formulation 
in terms of vector and scalar potentials. 


4.1. Starting Equations 

The following equations are used in this section: 
Continuity equation (Eq. 3.1.1) 


Dt 


pdiv v = 0 


4.1.1 


Cauchy* s equilibrium equation (Eq. 3.6.5) 

Dv 1 — 

— = f - grad li + Tgrad S + — div V 
Dt p 

Entropy evolution equation (Eq. 3.3.1) 

pT— = ^ ~ div q 
V Dt 

Also the expression for the acceleration given by Eq. 2. A. 4 


4.1.2 


4.1.3 


Dv dv 1 _ - 

— - — + -grad v - v x & 

Dt dt 2 


4.1.4 


will be used. 


It is understood that the equations of state, Eqs. 3.2.1, 
3.2.2 and 3.2.3 are available. It is also assumed that some 
constitutive relations (which give T and q in terms of other 
quantities) are available. Such relationships are discussed in 
Section 4A. The reason the introduction of the constitutive 
relations is postponed is to emphasize the generality of the 
vector/scalar-potential formulation, i.e., that the derivation of 
such a formulation is independent of the specific expressions 
used for T and q. 


4.2. Decomposition Theorem 

A fundamental theorem of vector field theory states that any 
vector field, in particular, in our case, the velocity field, can 
be decomposed as (see Serrin, p. 164-165) 

v = grad 4> + curl A 4.2.1 


4.1 



where <J> is called the scalar potential, whereas A satisfies the 
relation 

div A = 0 4.2.2 

and is called the vector potential. 

Note that the decomposition is not unique since any 
solenoidal irrotational field can be expressed either as grad <|> 
or as curl A . 

Taking the curl of Eq. 4.1.1 and using Eq. 4.2.1 and the 
vector formula 


curl(curl a) = grad div a — V*a 4.2.3 

one obtains that A must satisfy the Poisson equation 

V*A =-£ 4.2.4 

I _ 

Because of its relationship with the vorticity, curl A will be 
called vortical velocity 


v y = curl A 4.2.5 

Eq. 4.2.1 will be rewritten as 

V = grad 4> + V v 4.2.6 

In addition to the equation for ^an equation for (jj and one 

for 5 are required. These are derived in Sections 4.4 and 4.5 
respectively. Before doing that, however, an extension of 
Bernoulli's theorem is needed. In potential barotropic flows the 
equation for the potential is obtained by replacing Dp/Dt in the 
continuity equation with its expression obtained from Bernoulli's 
theorem. For viscous rotational flows such a theorem does not 
exist. Therefore a generalization of Bernoulli's theorem is 
presented in Section 4.3. 


4.3. Generalized Bernoullian Theorem 

The main innovation introduced in this work is the 
introduction of a generalized Bernoullian theorem. There exist 
several so-called Bernoullian theorems (Serrin, p. 153). The 
classical one is for unsteady irrotational inviscid barotropic 
(i.e., p=p(p)) flow in a conservative field (i.e., f=-gradJl , 

where fL is the potential energy). 

1^ / dp 

H = — + ~v-v +1 — + Q = H(t) 4.3.1 

at 2 Jp 


4.2 



Less known is the one for steady but rotational inviscid 
barotropic flow (Serrin, p. 153) 

+ 0 = constant over* Lamb surface 4.3.2 

(a Lamb surface is a surface defined by a network of vortex lines 
and streamlines). 


H 


- V. V + 
2 


f dp 

J P(P> 


It is apparent that Eq. 4.3.2 is considerably different from 
the classical Bernoulli’s theorem, Eq. 4.3.1, in that H is not 
constant in the whole field, but is a function of the location. 
Such a function is constant on the Lamb surfaces, or if one 
prefers, grad H is normal to the Lamb surfaces. The objective of 
this subsection is to obtain an expression similar to Eq. 4.3.2 
for viscous compressible flows, i.e., an expression which reduces 
to Eq. 4.3.2 for inviscid, barotropic steady flows. It is 
assumed in the rest of this work that -f is conservative so that 

f = -grad Q 4.3 .3 


In order to accomplish this, consider Cauchy's equation of 
motion, Eq. 4.1.2, which, using Eq. 4.1.4 and 4.3.3, may be 
rewritten as 

dv 1 

— + -grad v - v x f = 

at 2 


-grad ft - grad h + Tgrad S + -div V 4.3.4 

P 

or, using Eq. 4.2.6, 

04 > 1 

grad (— + — v - v + h + ft) = 
o t 2 


3 v. 


1 


” -r— -fvx£ + Tgrad S + -div V 
3t p 


4.3.5 


It is apparent that the right-hand side of Eq. 4.35 is 
irrotational (for instance by taking the curl of both sides). 
Hence the integral of minus the (dissipative) terms on the right 
hand side. 


f ( 

T> 


-J 


* 3v. 


(-— 1 - v x 5- Tgrad S - -div V )• 4.3.6 

at p 

is path independent (suitable branch surfaces are introduced for 
multiply connected regions), i.e., fp is only a function of x and 
is such that 


3v v _ - 1 

grad f = v x ? - Tgrad S - -div V 

o dt p 


4.3.7 


Combining Eqs. 4.3.5 and 4.3.7 yields 


3<p 1_ _ 

— + — v • v + h + 
3t 2 


A + f 0 = Mt) 


4.3.8 


4.3 



where n o is a function of time (for instance, if the fluid at x 
= -® is in uniform translation with speed U m , IT = U 1 /2 + h + 
fp M where h^ and fp^are the values of h and fp at Xj = -»). 
Equation 4.3.8 may be written as: 

H= + A v 2 + h +il = constant (over fp-surface) 4.3.9 
d t ^ 

(where a f^-surface is a surface defined by f^ — constant). 

Equation 4.39 is the generalized Bernoullian theorem. It is 
apparent that it reduces to Eq. 4.3.2 for steady, isentropic 
(S=constant), inviscid (V = 0) flows, since, in this case, 

H=H 0 -f D 4.3.10 

is constant on the Lamb surfaces (for, vx£is by definition 
normal to Lamb surfaces). 

It is understood that Eq. 4.3.8 is a 'formal' (rather than a 
substantial) generalization of Bernoull i's theorem. In other 
words, no claim is made here a^out the discovery of a new 
physical concept. More appropriately, Eq. 4.3.6 should be 
thought of as a convenient formal expression that can be used to 
obtain a differential equation for the scalar potential (see 
Section 4.4). The determination of the function fp(x) is a 
problem as complicated as the original one, unless the hypotheses 
under which Eq. 4.3.2 (isentropic inviscid steady flow) is valid 
are satisfied. However, it will be seen in Section 4.4 that, at 
least for steady state, the evaluation of f^ i s not necessary 
(see Eq. 4.4.10). Therefore, at least in that sense, the above 
generalization is a powerful tool because it allows for the 
derivation of the differential equation for the scalar potential 
4>. Also, it may be worth noting that, integrating along a 
streamline. 

x 

f n (x) = ( (— ’l - T grad S -- div v)*dx 4.3.11 
J db P J 


4.4. Differential Equation for Potential 


Consider the continuity equation, Eq. 4.1.1, which may be 
rewritten as 

1 Dp 

div v = - — — 4.4.1 

p Dt 

Combining with Eqs. 4.2.1 and 4.2.2 yields 

2 i 1 Dp 

V*4> = - - -1 4.4.2 

p Dt 

Assuming h and S as fundamental variables of state, one obtains 

.dp, Dh dp, 


1 


DS 



4.4.4 


Note that (see Eq. 4.C.3) 


A d p | _ 3p. _ 1 
p ah *s 3p'S 

where a is the usual isentropic speed of sound* 


Introducing for notational convenience 

P »s h 


Eq. 4.4.3 may be rewritten as 


7*+ = - i - . 


DS 


2 Dt 


Dt 


4.4.5 


4.4.6 


The enthalpy h can be obtained from the generalized 
Bernoullian theorem, Eq. 4.3.8, to yield 

1 D 34> 1 _ . v DS 

" 9 = — ( — + -v-v + f + n) + B— 

a *Dt at 2 » Dt 


1 D „ DS 

= ( — + v- v + f + Q) + B 


4.4.7 


a 2 Dt 3t Dt 

(where c indicates that v is kept constant while applying the 
operator D/Dt), or 

„ 2 . 1 » C + _ Df B j>n DS 

7 4> = — ( + v + — + — )+ B — 4.4.8 

a* Dt 2 Dt Dt Dt 

where in D* /Dt 2 the term v in the first substantial derivative 
is kept constant daring the second material differentiation. 


Note that, according to Eq. 4.3.7 
Df a _ flf 3v v 

Bt " dt + v ‘ grad f D" + v *Tt ” Tf ' grad S ~j v>divV 4 * 4 ’ 9 

which indicates that, as mentioned above, in steady state the 
explicit evaluation of f D is not neC essary. 


4.5. Vorticitv Dynamics Equation 

In order to complete the formulation an equation for is 
needed. Taking the curl of Eq. 4.3.4 yields 

dv i _ 

curl— - curl(v x?) = grad T x grad S + curl(-divV ) 4.5.1 

ot p 

Noting that 

curl(a x b) =-a.grad b + b.grad a + adiv b - bdiv a 4.5.2 


one obtains 


4.5 



cur 1 ( v x £ ) =-v- grad 5 + S'* grad v + vdiv£ - 3 div v 


— — c 1 Dp “T 

= -v • ^r ad 5 4 b ' grad v 4 ” ^ 

and, combining with Eq. 4.5.1 ana noting- that 

d v 6 dZ 

curl--=-- curl v= ~ 
dt dt dt 


one obtains 


sz x - 

— + v- 
dt 


— — Dp y 

grad 5 - grad v — S 


p Dt 


4.5.3 


4.5.4 


= grad T x grad S + curl(-div V) 

P 


4.5.5 


or 


— (?) = — . grad v + -grad T x grad S + -curl(-div V) 4.5.6 
Dt p p p p p 


which is an extension to unsteady viscous fluid of Vazsonyi 
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vorticity dynamics equation (Serrin , p. 189). 
4.6. Summary of results 


The results can be summarized as follows. The fundamental 
unknowns are the potential <f> , the vorticity 5, and the entropy S. 
The corresponding equations are: 


Potential (Eq. 4.4.8) 

~ 1 D c 2 £ _ D % Df D DS 

V — — ( + v • - — — + — — + — — )+ D — 

a 2 Dt* Dt Dt Dt Dt 


4.6.1 


Vor t i c i.ty (Eq. 4.5.6) 


D 


-(-) = -.grad v + -grad T x grad S + -curl(-div V) 4.6.2 
Dt p p p * P P 


Entropy (Eq. 4.1.3) 

DS 

pT— = $ - div q 

Dt H 


4.6.3 


In addition the enthalpy is given by the generalized 
Bernoullian theorem (Eq. 4.3.8) 


d4> 1 

— + — v . v + h + IT + f n = H 0 
6t 2 D 


4.6.4 


whereas the velocity is given by Eq. 4.2.1 



v = grad<£ + curl A 
where A satisfies Eq. 4.2.4 


4.6.5 


V*A 


4.6.6 


The above equations have more unknowns than equations and 
require equations of state giving p,p and T as functions of h and 
S (such as those for an ideal gas considered in Section 4. A) and 
constitutive relations j[or the heat flux vector, q and the 
viscous stress tensor, V, such as the Fourier law for heat 
conduction and the Cauchy-Foi sson law of viscosity (which yields 
the well known Nav i e r-S t oke s equations of motion) which are 
discussed in Section 4.E. 

The above equations with appropriate state equations and 
constitutive equations may be solved approximately as discussed 
in Section 5 for both attached and separated flows. The possible 
advantages of the scalar/vector potential formulation over the 
primitive - variable solution of Navier-Stokes equation are also 
discussed in Section 5. 


It may be noted that for inviscid (V = 0), adiabatic (q = 0) 
flows, in the absence of external forces (fl= 0), the above 
equations reduce to 




a 



+ 




4,6.7 


D ? ? _ 
— (-) = -.grad v 
Dt p P 


1 

- grad T x grad S 
P 


4.6.8 


DS 

Dt 


= 0 


4,6.9 


Note that integrating Eq. 4.3.6 along a streamline (using dx = 
(v/lv|)ds = v d9, where s is the arclength along a streamline 
and de = ds/ I vl ) 


-I 


3 v 

. t £ 5 , d9 
at at 


4.6.10 


since 


DS 

Dt 


as 

— — + v .grad S = 0 


4.6.11 


In particular, for steady state, noting that (see Eq. 4.3.7) 
v.grad f D = 0, 


ijtj 1 - , v 

v <p = — v. grad - 

2 2 . 

a A 

— * «J^ 

v - grad(— ) = —.grad v + — grad T x grad S 
P P P 


4.6.12 


4.6.13 


4.7 



4.6.14 


-v. v+h = constant on Lamb surfaces 
2 

S = constant along streamlines 

4 .A. Equat ion o f State for an Ideal Gas 

In this section the necessary equations of state for an 
ideal gas are presented. An ideal gas is a gas which satisfies 
the equation 

- = P t = RT 4.A.1 

P 

where R is a constant. 


This implies (see Eq. 3. A. 6) 
de * 

Or 


= 0 


and therefore 


r 

y 


e = e (T) = f c v (T) dT 


where 


de 

° v ” dT 

In addition (see Eq. 3.6.4 and 4.A.1) 

T 

h = li(T) = J c p (T) dT 

where 


c = c + R 
P v 

Finally (see Eq. 3.2.4) 
de p 

dS = — + — dx 
T T 


or 


; - s ° ■ 1 


dT 

c — + R lnx 

V f 


4. A. 2 


4. A. 3 


4. A. 4 


4. A. 5 


4. A. 6 


4. A. 7 


4. A. 8 


Equation 4. A. 5 can be solved for T to yield T = T(h). This can be 
used in Eq. 4. A. 8 to yield z = x(S,h) and finally Eq. 4.A.1 can 
be used to obtain p = pRT = p(S,h) which are the three desired 
equations of state relating temperature, density arid pressure to 
the fundamental state variables S and h. 

In particular, for ideal gases with constant coefficients, 

i. e . , 


4.8 



constant 


4. A. 9 


c v " 

and hence (see Eq. 4. A. 6) 

Cp ■» constant 

Eqs. 4. A. 3 and 4. A. 5 yield 
e = c v T 
h = c p T 

(where the constant of integration has been avoided by 
e = h = 0 for T = 0) , whereas Eq. 4. A. 7 yields 

dT dr dp dr 

dS = c v -- + r* — = c — + c„ — 


p p X 


or 


pr*= e S/c - 

(where the constant of integration has been avoided by 
S = 0 for p = r = 1) and 

RT t*" 1 = e 8 / 0 * 

Hence, one obtains explicit expressions 
T = T(h) = h/c 

P 

t = r ( h, S) = (R he Sfc */c = 

P 

* 

p =p(h, S) = Ch lli) zr/(y-iJ e -s/B 


It may be worth noting that (see Eq. 4.4.5) 

. i?J| „ i 

p 3Slh x 0S|h R 


4.A.10 

4. A. 11 
4. A. 12 
choosing* 

4. A. 13 

4. A. 14 
choosing* 

4. A. 15 

4. A. 16 

4. A. 17 
4. A. 18 

4. A. 19 


is a constant. Also, note that the fundamental equation of 
state, Eq. 3.2.1, is (see Eq. 4. A. 12) 


e 


- i e s/c v 

ar-i ri - 1 


4. A. 20 


i.e., the formulation for ideal gases with constant coefficients 
is equivalent to postulating the above equation for the 
fundamental equation of state. 

* This is consistent with the third law of thermodynamics. 


4.9 



4 . D. Const itut ive Relations 


In this section, the constitutive relations are discussed: 
in particular it is indicated how the Cauchy-Pois son law of 
viscosity (see Eq. 3.3.4) 

V = X0I +3.^5 4.B.1 

(where 0 = «<*v v, and the viscosity coefficients X and jx. are 

positive functions of thermodynamic variables) and the Fourier law 
of heat conduction 

q = -k grad T 4.D.2 

(where the conductivity coefficient k is a positive function of 
the thermodynamic variables) may be obtained by introducing 
certain additional restrictive (but 'reasonable') postulates. 

The first postulate introduced here is that the entropy 
condition, Eq. 3.4.4, be satisfied by each individual term (see 
Eq . 3.3.6) 

$ = V : D > O 4.E.3 


and 


q.^rad T > O 4.B.4 

These are the mathematical expressions for the familiar 
statements that mechanical deformation dissipates energy 
(transforming it into heat), and that the heat flows in the 
direction of the temperature gradient. 

Next, consider a set of four postulates introduced by Stokes 
(see SerrinJ p. 231) 

1. V is a continuous function of the deformation tensor D 
and is independent of all other kinematic quantities. 

2. V does not depend explicitly on the position x (spacial 
homogeneity) . 

3. There is no preferred direction in space (isotropy). 

4. When D = 0, V = 0 as well. 

A fluid satisfying the above postulates is called Stokesian. 
Mathematically speaking, the first, second and fourth postulates 
imply 


V = f(D) ( f ( 0 ) =0 ) 4.D.5 

(The dependence upon two thermodynamic variables, such as p and T, 
while not stated explicitly, is not being excluded). The third 
postulate implies that Eq. 4.B.3 is invariant under all 
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rectangular coordinate transformations. 

In addition, consider another postulate: 
5. The function f(D)is linear in V. 


These postulates are sufficient to obtain the Cauchy- 
Poisson law of^viscosity, Eq. 4.B.I. (A proof of this equation is 
given in Serrin T pp. 233-234). It nay be worth noting that if the 
fifth postulate, which is the most restrictive one, is removed 
then, as shown in Serrin pp. 231-232, f(D) must be of the type 

T=«<I + pD + YP 2 4.B.6 

where a, p and y are functions of the invariants of D (as well 
as of the thermodynamic variables). 

Introducing similar postulates for q yields the Fourier law 
for heat conduction, Eq. 4.E.2. 

4 . C . Speed of Sound 


Note that 


dp 


= i£| 

3hls 


dh + 


dpi 

0S(h 


dS 


4.C.1 


On the other hand, using Eq. 3,6.2 


dp 



^(pdh-pTdSW^ 


dS 


dp I . dp . dp , 

= P — dh + (- — pT + )dS 
®P*S °P's *P 


4.C.2 


Comparing Eqs. 4.C.1 and 4.C.2 yields 

- I _ £P| _ i 

p 3h] 3p l a a. 


4.C.3 


where a is the usual (isentropic) speed of sound, defined as 

4 = ‘fplj’ 4 - c - 4 
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S ect ion 5^ 

Cone luding Remarks 


5.1. Comments 


A general formulation for viscous, compressible flows has 
been presented. In order to discuss the advantages of the 
formulation, consider the case of the attached flow around an 
isolated wing at a high Reynolds number. In this case the region 
of nonzero vorticity is limited to a thin region around the 
wing (boundary layer) and a region behind the wing (wake), and 
the numerical results for the pressure distribution obtained 
under the assumption of isentropic irrotational flow are 
generally in excellent agreement with the experimental ones, both 
in the subsonic and supersonic regimes. This would indicate that 
the effects of the presence of the vorticity in the field are 
small and therefore the solution to Navier-Stokes equations may 
be obtained as a correction of the potential flow formulation. 
In order to better understand this, consider the relationship 
between the inviscid adiabatic formulation (Eqs. 4.6.7 to 
4.6.11) and the potential formulation. If the flow is initially 
isentropic (e.g., at rest at t = 0), Eq. 4.6.9 yields that S = 
S Q (where S Q is a constant) at all times. Then if the flow is 
initially irrotational (e.g., at rest at t = 0), Eq. 4.6.8 yields 
that 5= 0 at all times. Therefore v y = 0 at all times, and Eq. 
4.6.10 yields fp = 0 at all times, and Eq. 4.6.7 reduces to the 
classical formulation for the velocity potential. 

Next consider the relationship between inviscid adiabatic 
formulation (Eqs. 4.6.7 to 4.6.11) and the formulation for 
viscous conductive flows (Eqs. 4.6.1 to 4.6.6). In both cases 
we will assume that the flow is initially at rest. If the 
viscosity and conductivity coefficients are snail (i.e., high 
Reynolds number flows), then Eq. 4.6.3 yields that S is 
approximately constant in the outer region (the fluid volume 
minus the boundary layer and the wake region). In addition Eq. 
4.6.2 yields that the vorticity is approximately equal to zero in 
the outer region. Next consider the velocity, v v , induced by 
the vorticity (see Eq. 4.6.6). With a suitable choice for the 
boundary conditions for Eq. 4.2.4, it is possible to obtain that 
A = 0 (approximately) in the outer region, since the velocity 
obtained from Eq. 4.6.6 is irrotational in the outer region. As 
a consequence of the above remarks, 3S/ 3t and Bv / 9t are zero 
in the outer region and so is fp. Therefore, the additional 
(nonpotential) terms in Eq. 4.6.1 are all approximately equal to 
zero in the outer region and hence the only modification in the 
equation for the velocity potential due to the effect of 
viscosity is the presence of some apparent sources of mass in the 
region of the boundary layer and the wake. In this case an 
integral equation formulation would be ideal to approach this 
problem (at least in the subsonic and supersonic regimes). 
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5.2. Recommendation for future work 

A general formulation for scalar/vector-potential 
decoiry^sit ion for viscous compressible flows has been presented. 
This formulation should not be considered the last word on the 
subject but merely the first. Several issues have not been 
addressed here, such as separated flows, turbulence, and the 
presence of shock waves in the field and the vorticity generated 
by them. However before addressing such issues it seems 
appropriate to assess the formulation presented here by 
developing a numerical algorithm (and corresponding computer 
program) and comparing the results against existing numerical and 
experimental ones. 
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